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Abstract 

We investigate cosmological models with a free scalar field and 
a viscous fluid. We find exact solutions for a linear and nonlinear 
viscosity pressure. Both yield singular and bouncing solutions. In the 
first regime, a de Sitter stage is asymptotically stable, while in the 
second case we find power-law evolutions for vanishing cosmological 
constant. 
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1 Introduction 

Studies of the dynamics of a scalar field together with other forms of matter 
usually assume that this behaves as a perfect fluid. Few authors have taken 
into account dissipative effects [|TJ [||] j|, even though processes like par- 
ticle creation might have been important in the early universe. In this paper 
we study the evolution of a universe filled with a massless minimally coupled 
scalar field and a viscous fluid. We find exact solutions of the Einstein equa- 
tions in a Robert son- Walker metric and we analyse their asymptotic stability 
by means of the Lyapunov method ||. 



2 The model 

We investigate the evolution of a universe filed with a scalar field and a 
viscous fluid. The scalar field is free and minimally coupled to gravity, so 
that it obeys the equation 0<f> = 0. In the case of the homogeneous isotropic 
Robert son- Walker metric 



ds 2 = dt 2 - a 2 (t) 



dr 2 



kr'' 



+ r 2 (d9 2 + sin : 



(1) 



the scalar field equation becomes 



4> + 3H<j) = (2) 

Besides, only the bulk viscosity needs to be considered. Thus we replace 
in the Einstein equations the equilibrium pressure p by an effective pressure 



3H 2 



2H + 3H' 



^ 2 + p-3^ + A 
2 er 
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where H = a/a, ' = d/dt, p is the energy density, a is the viscous pressure, 
A is the cosmo logical constant and we use units c = 8irG = 1. As equation 
of state we take 



P = (7 - !)P 
2 



(5) 



with a constant adiabatic index < 7 < 2, and we assume that a has the 
constitutive equation 

a = -3(H (6) 

where ( > is the bulk viscosity coefficient. Thus, we must solve equation 
(H) together with the Einstein equations Equation (§) has the first 

integral 

<P = ~s (7) 

cr 

where C is an arbitrary integration constant, and we are interested in the 
case that C 7^ 0. Then, adding (Q) and (f|), and eliminating p, we arrive at 

2H + 3 7 # 2 - 3(H = + (2 - 3 7 ) - 2 + 7 A (8) 

To find exact solutions of this equation we make the change of variable a = 
v u , and we consider two cases: ( constant and ( ~ H. We assume now that 
C is a constant. This is a good approximation when the thermodynamical 
variables do not change too much. This has been studied in several papers 
[0] (1 i Q- Choosing v = 2/(37), the equation (§) becomes 

v + ^(v + Mv m + Nv n --^-v = (9) 

where M = (3/8) 7 (2 - 7 ) C 2 , m = 1 - 4/ 7 , N = (3/4) 7 (3 7 - 2)k, n = 
1 — 4/(37). ^ is linear for 7 = 2 and k = 0, so that we can find its general 
solution. In this case, the evolution of the scale factor is the same as in 
a model without the scalar field and a bulk viscosity coefficient (/2 [0 [§]. 
Thus, we just quote the solutions without further analysis. 

Let us consider first that A = 0. We find 



a(t) = 



Aexp (Ja) 



2C 

A0M = ^ln 



exp 



v4exp(|Ct 



nl/3 

B 

m 



-B 



(10) 
(11) 
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/N |C 2 A 2 exp(3a) -C 2 , N 

2[Aexp (fa) + 5 

where A and -B are arbitrary integration constants. The requirement p > 
implies that there is a minimum time for the validity of this solution. Thus it 
describes a geodesically incomplete manifold, which in this sense is singular 
E|. Clearly, in its asymptotically de Sitter regime it is physically justified 



the approximation of constant £. However, for shorter times, this assumption 
may not be justified. 

The asymptotically de Sitter stage occurs for any other value of 7 and 
k. In effect, assuming a ~ exp(Hot) for t — > 00, with i?o a constant, we find 
that Hq = C/7- To verify the asymptotic stability of this solution we use a as 
the independent variable and we take the Liapunov function L = (H — Hq) 2 . 
Then, we find that 

L' = —37— + O (^-) (13) 



a \a 3 



is negative for large times. 

Let us consider now that A ^ 0. There is a critical value for the cosmo- 
logical constant A = — 3( 2 /16. Thus, we must distinguish several cases: 
Two kinds of solutions appear for A > Ao- Singular evolutions 

a(t) = a e 3A ' /4 [sinh (uoAt)} 1/3 (14) 
or nonsingular evolutions 

a(t) = a e 3At/4 [cosh (ooAt)} 1/3 (15) 

where u = [3 (A — Ao)] 1//2 , At = t — to and ao, to are integration constants. 
For the scalar field we find 



While p — > (3/8)C 2 + (3/2)(u > for t — > 00, in the nonsingular solutions 
p becomes negative before a minimum time. On the other hand, for singular 
solutions, we find 
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Ai -° (17) 

so that the energy density may remain positive along all the evolution. 
In the case A = Ao, the evolution is always singular 

a(t) = a (At) 1 / 3 e 3A */4 ( 18 ) 
In the case A < A , the scale factor recollapses to a second singularity 

a(t) = a e 3Ai/4 [sin (\oj\ Ai)] 1/3 (19) 

All singular evolutions have particle horizons because a ~ At 1/3 as At -> 0. 

For other values of 7 and k, there is also a critical value of the cosmological 
constant A = — 3£ 2 /(47 2 ), so that for A > A an asymptotically de Sitter 
evolution occurs for two values of H 



2 7 



4 \ 1,2 

C ± ( C + ~7^ 



(20) 



We verify that the behavior a ~ exp(H + t) is asymptotically stable by means 
of the Liapunov function L + = (H — H + ) 2 , which satisfies L' + < for large 
times when H > 0. 

Nonlinear viscous effects has been shown to arise as a phenomenological 
description of the effect of particle creation in the early universe []13| |, and 



cosmological models with this kind of fluids has been considered in JTj] and 

Here, we assume that ( = aH, with a constant a such that 7 > a > 0, 
and we consider only expanding evolutions. Following the same procedure 
as before, we arrive at the equation 

v + Mv m + Nv n - ^-v = (21) 

where now v = 2/[3(7 - a)}, M = (3/8)(7 - a) (2 - 7), m = 1 - 4/( 7 - a), 
N = (3/4)(7 — a) (37 — 2)k and n — 1 — 4/[3(7 — a)]. As this case behaves 
as a conservative mechanical system, we may obtain (at least formally) its 
general solution: 
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At = I , dv (22) 



'2 [E - V(v 
where E is an integration constant and 

V{V ) = -4-,-" + 4-„»« - %J (23) 
m + 1 n + 1 4z/ 3 

V(v) = — — — t> m+1 + AHnf — ^-^t> 2 7 -a = - (24) 
v ; m + 1 4i/ 3 v ; 



A qualitative analysis of (|22| ) shows that expanding evolutions may begin 
either at a singularity or a bounce. Both bounded and unbounded singular 
solutions occur, but only for A < are there bouncing solutions that reach 
a maximum. 

We quote two cases for which the equation fl2~T|) becomes linear, so that we 
may obtain its general solution in closed form (for simplicity we take A = 
). The first one is 7 = 2 and k = 

a{t) = (AAt) u (25) 

AMt) = -^p (26) 

vv ; A 3v (l - 3i/) V ; 

p{t) = i^~ 2(iMr (27) 

The other solution is for 7 = 2, a = 2/3 and k 7^ 0. We find 

a(t) = \AAt - 2kAt 2 \ 1/2 (28) 
* * / \ 2C A - 4k At 

A<P{t) = - — m 29 

A 2 \AAt - 2kAt 2 \ 1/2 

3 (A - 4kAt) 2 1 C 2 3fc 

PW ~ 4\AAt-2kAt 2 \ 2 2\AAt-2kAt 2 \ 3+ \AAt-2kAt 2 \ (,>) 

We find also solutions of the form a = aoAf for special values of a: a. 7 = 
2, a = 1, a = (4/3)(l+£;/a 2 ); b. k = 0, a = 1/3, a = 7+3(2-a)C 2 /(2a§)-2. 
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3 Conclusions 



We have found exact solutions of the Einstein equations with a free scalar 
field and a viscous fluid source, in a homogeneous isotropic metric. We have 
considered cosmological models with two forms of the bulk viscosity pressure: 
linear, with a constant coefficient, and cuadratic in the Hubble variable. For 
the first case the picture is similar to the case without scalar field and we 
show that the de Sitter evolution is asymptotically stable. Also, we find 
that the physical requirement of positivity of the energy density of the fluid 
makes some of these solutions geodesically incomplete. For the case of a 
nonlinear viscous pressure we reduce the equations of this model to that 
of a conservative mechanical system. Thus, we are able to give its general 
solution in implicit form, and besides we show several cases for which an 
explicit solution arise. Both, singular and bouncing solutions arise, and we 
find cases in which the scale factor reach a maximum or grow without bound. 
In the latter case, when A = 0, the scale factor has a power law behavior for 
large times. 

Dissipative effects like particle production may have been important in 
the early universe, for instance in the reheating period at the end of the infla- 
tionary era. It is generally assumed that this stage of accelerated evolution 
was driven by a self-interacting scalar field. Thus, we consider of interest 
to investigate further the models of this paper, taking into account also an 
interaction potential. 
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